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Abstract 

The solution of the classical open-chain n-body Toda problem is derived 
from an ansatz and is found to have a highly symmetric form. The proof 
requires an unusual identity involving Vandermonde determinants. The ex¬ 
plicit transformation to action-angle variables is exhibited. 
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The Toda chain is one of the paradigmatic examples of an integrable many-body 
system of interacting particles. The discovery of its conserved integrals of motion]^, 
1^ and its snbseqnent solntion[P, ^ were important steps in the development 
of the theory of integrable systems|P]. An almost universal feature of analytical 
studies of the Toda system is the use of the Lax pair formalism. In this paper, 
an alternative derivation of the solution of the classical open-chain n-body Toda 
system is given. 

The derivation proceeds essentially from an ansatz about the form of the solu¬ 
tion and therefore lacks the power and generality of the Lax pair treatment. The 
solution however has an elegant structure which is not evident in previous repre¬ 
sentations. More, it can be interpreted as the classical canonical transformation 
from the Toda system to a free theory. This is an important clue to construct¬ 
ing the classical and quantum solutions by a sequence of elementary canonical 
transformations!^. Following the successful solution of the three-body Toda prob¬ 
lem with this approach!^, work is in progress on the classical and quantum open- 
chain n-body problems. 

The Hamiltonian for the (n-l-l)-body open chain Toda system is 

n+1 n 

H = + ( 1 ) 

k=l k=l 

The arguments of the exponential potentials can be interpreted as expressions for 
the root vectors of An in the Cartan basis^]. A coordinate transformation will 
put the root vectors into the Chevalley basis and separate out the motion of the 
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center of mass. The transformation is given by 


?! 


1 Qn+l 

qi+ , 
n + 1 

Qk 


Qk-i + qk+ , {2<k<n) 

n + 1 

Qn+1 


1 9n+l 

9n+ , 

n + 1 




Pk 


„ , 1 ( E (^ + 1 j)Pj E ) + Pn+l 

j=l ^ 

Pn+1 




The transformed Hamiltonian is 


( 2 ) 


(1 < A: < n). 


n k—1 


+ k)pi + Y.Y1 2 j (^ + 1 - k)pjPk ] + 

Z{nl) \ A:=2i=l 


n—1 


+e 


2qi-q2 


^qk-qk-l-qk+1 _j_ 


(3) 


k=2 


This leads to the equations of motion 


Qk 

Qn 


_^2qi-q2 
_Qfc — 1 ^fc + l 

_ Qn—l 


{2 <k <n — l) 


(4) 


The solution of these equations has the remarkably simple form 

n+1 

where A^(ji, • • • ,im) is the square of the Vandermonde determinant 

A (jl, . . . ,jm) = {pji ~ hj'fc) ) (6) 

ji 
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and fk and fik are arbitrary constants, satisfying 


n+l 

H/a: = A-2(l,...,n + l), (7) 

k=l 

n+l 

= 0 . 

k=l 

(There are additional constraints on the range of the + if one requires the be 
real.) The solution has 2n free parameters as required. The solution in the original 
variables is determined from the transformation (§) to be composed of ratios of 
these solutions times a factor for the center of mass motion. 

To derive the solution, make the ansatz 


e 


n+l 

E f. . 


( 8 ) 


where the are arbitrary real numbers. Note that this ansatz defines a variable 


g-<?n+i ^ 


(9) 


Such a variable might naturally appear in the hnal equation of to give + = 

variable does not appear can be interpreted as meaning 
that one has set g„+i = 0. This ultimately is the origin of the restrictions (|^) on 
the fk and /i^. The equation for gi is also of the form of the others if there is a 
go = 0. The open-chain Toda system thus has hxed endpoints in this sense. The 
ansatz and solution are compatible with the slightly more general problem where 
g-g„+i ^ Then, in the solution, one would have Ilfcii fk = cA“^(l,..., n -1- 1) 
and Yfkti l^k = «• 
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Consider e Differentiating twice and multiplying by e leads to 


gmC = e - ((9te 


(10) 


But from the equations of motion — g^e _ g qm-i qm+i (fusing ^iq = 0 = g^+i). 
Substituting the ansatz into the resulting equation gives (2 < m < n — 1) 


n+l 




(11) 


k^< -<km 
Jl<fcl 


i=l 


i=l 


n+l 


= V f- ■ h k 


3l<---<3m-l 

fcq < •• •<fc77T,_)_l 


The equation for m = 1 is 


n+l n+l 

E fhMPh - +'■»>' = E (12) 

ii<i2 ii<t2 

The equation for m = n involves where 1 < ji <■■■<+< ''^ + T As one 

is choosing n integers out of n + 1, this is more succinctly labelled by ff where r 
is the integer which is not in the set. Similarly means the two integers r ^ s 
do not appear, and the indices of / are the remaining n — 1 integers. With this 
notation, the equation for m = n is 


^ _i_1 

E fM-k. + «)< 

r<s 


^ f^Qi-^kr-^ks+J2lll^kk)t_ 
r<s 


(13) 


Assume the are all distinct and that they have no accidental degeneracies in 
their linear combinations, such as fij^ + + fij^. The asymptotic behavior 

of the exponentials can be used to equate like terms in the sums. The degenerate 
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cases can be recovered later by continuity in the Hk- Let 



(14) 


where the fj^, are (so far) arbitrary constants and A^(ji,... ,jm) is the square of 


the Vandermonde determinant (P). With this dehnition, the m = 1 equation (|T^ 


is easily verihed. The m = n equation (|T^ is satished if the constraints (|^ on the 
fk and /ifc are imposed. The proof that Eq. (|TT]) is satished reduces to a heirarchy 
of identities for Vandermonde determinants. 

On the left hand side of (0), there are two sets of indices {ja} and {k^}. 
Since ji < ki, at most they can have m — 1 indices in common. The asymptotic 
behavior of the exponential is given by a sum over the /i* indexed by the combined 
set S = {ia^kfj}. Different partitions of S into sets {ja} and {kjs} {ji < ki) will 
have the same asymptotic behavior. The number of such terms will depend on 
the number of distinct indices between the two sets, and these constitute separate 
cases. Let 2r denote the number of distinct indices. 

Consider the case r = 1, labelling the common indices si,..., and the 
distinct ones ji and k (ji < k). There is a unique term on both sides of (|TTD with 
the asymptotic behavior given by this set of indices, and one has 



(15) 




where the curly brackets indicate the indices should be arranged in increasing 


order. Using (p!4D, the constant factors /j cancel and one has a relation between 
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Vandermonde determinants 


^ (jl; '^1? • • • 7 l)^ (^7 -^17 • • • 7 


(16) 


A (^i, • • • , 577 j,_i)A (ji, Sij • • • 7 1)‘ 


Using the relation 

m—1 

A^(A:, si,..., Sm_i) = A 2 (si, ..., s^_i) Y[ (/^fc - (17) 

and its relatives, the dependence on the common indices is seen to cancel and one 
is left with the identity 

f^kf = (18) 

It is a general feature for all r that the dependence on the constant factors and 
the common indices cancels on both sides, so without loss of generality one can 
focus on the distinct indices alone. Reindex the set S of distinct indices by the 
integers 1 to 2r. Partition S into two sets a = {1, Q! 2 , • • •, ctr} and [3 = {/3i,..., Pr} 
and denote the collection of such partitions Pap. Separately partition S into sets 
7 = { 71 ,..., 7 r-i} and 5 = {(5i,..., calling the collection of partitions P.yp. 

Denote A^(Q;;r) = A^(l, 02 , • • •, a^) and similarly for the rest. The number r 
of indices involved in the Vandermonde determinant is made explicit to reduce 
confusion. Both sides of Eq. (O) will be equal if the following identity between 
Vandermonde determinants holds 


^ A2(a; r)A2(/3; r)(^ = Y. r - l)A2(h; r + 1). (19) 

PaP “ /3 P^5 
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It seems likely that this identity has a group theoretical interpretation, but in 
its absence, the identity can be proved inductively as follows^. Divide both sides 
by 2r). This gives the equation 



( 20 ) 


Denote the left hand side of the equation by L^. and the right hand side by Rr- 
The equation Li = Ri holds trivially. Assume that L^-i = Rr-i- The inductive 
step will be made by considering the pole structure of Lr and Rr. Since Lr and 
Rr are analytic functions of the /Uj without zeroes, if they can be shown to have 
the same residue at all of their poles, they must be equal. 

Choose two indices from the set S, neither equal to 1, and let their associated 
fii be labelled 2 : and a. (The index 1 is special because it has a preferred role in 
the partitioning. Which of the original /ij is associated to the index 1 is however 
arbitrary, so one can investigate the pole structure at the /Xj missed here by rein¬ 
dexing the set S.) Let S' denote the set S with these two indices removed, and let 
a', (3', 7 ', 6' denote partitions of S' as dehned above with r replaced by r — 1 . 

Consider the residue of z = a. Lr has a double pole at z = a if z ^ a 

and a & j3 OT vice versa. In the former case, the residue is computed to be 


Res L 


z=a 



'r 


z£a,a£l3 


X/a' La' ■ — ^ ^ p. , ^ . I- / 



8 









In the alternative case z ^ P, a ^ a, the residue is 


Res Lj 

z=a 


z£l3,a£a 

2 


E 


E 


(22) 

1 \ (Sa' ho' S/3' h/3') 


Os'(fl Pa'/3' ^ ^a' f^a' S/3'h/3' a' ® t^a'^ Y\a' ,13'if^a' I^f3')'^ 

Adding these, the residue of L,. at z = a is 

2,L^ — 1 


ResLr = 


-E—=—• 

ns'(“ - ^ a- li-S' 


(23) 


The residue at 2 ; = a of R,. is similarly composed of terms where z G a E S 
and vice versa. The full residue is 


ResRr = 


2Rr_i 1 


(24) 


ns'(“ - y a - /Is/ 

This is seen to equal the residue of L^. a.t z = a, given L^-i = Rr-i- Since this 
result holds for all pairs of the original /Xj, one concludes that = Rr and the 
induction is complete. 

To exhibit the solution as a canonical transformation from (^) to a Hamil¬ 
tonian independent of coordinates, one must introduce hnal coordinates and mo¬ 
menta and hnd a relation between them and the fj and fik so that the transfor¬ 
mation is canonical. It is clear that one can redehne fj by an overall constant. 




The arguments of the exponentials then dehne the hnal coordinates 


(25) 


Xj — fJjjt Xj. 


(26) 
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There should only be n independent degrees of freedom, and the coordinate Xn+i = 
— Yh=i Xi is not independent because it is related to the others by the constraints 
(^. It is useful however to introduce a temporary form of the final Hamiltonian 


1 n+l 

= 2 E f*?' (27) 

1=1 

The Pj are not the momenta conjugate to Xj because if the the constraint //„+! = 
— IJ'i were eliminated, the wrong Xj would follow from Hamilton’s equations. 
It is necessary to introduce n momenta kj conjugate to the Xj, so that ij = = 

Hj. The relation between kj and pj is found to be 

n 

kj fJjj /ij ^ ^ fJ^i ( 28 ) 

i=l 

or in reverse {j ^ n + l) 

1 

~ “ I 1 X! 

“I" i=i 
1 " 

P-n+l = I ~ ki 

n+1^ 

The final Hamiltonian is then 


r? ” 1 

2{n+l)+, ^ n+l + ^ ^ 


(30) 


The next step is to find an equation for the evolution of the original momenta. 
This is easily done by taking a time derivative of the solution 


n+l 




(31) 
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to find 


n+1 _ 

- = Y. fh--- • • •, Jm)(/^ji + ■ ■ ■ + (32) 

Using Hamilton’s eqnations with the Hamiltonian (^, one can express Qm in terms 
of the momenta as 


m—1 


Qm = + 1 - m)pm + Y ^ ~ ^)Pi + nr(n + 1(33) 

yTl -|- j i=m+l 


The result is 


m —1 n 

m{n + 1 — m)pm + Y^ + 1 — rn)pi + Y, nT,(n + 1 — z)pi]e“'^'" = 

2=1 2=m+l 

n +1 

= X! hi'" ■ • • ljra){Pjl H-h • (34) 

Finally, by reqniring that the Poisson brackets be preserved nnder the transfor¬ 
mation, one can determine the fj in terms of the kiS. The result is that {j ^ n-|-l) 

n 

fj = {-iy-^k-^Y[{kj - ki)-\ (35) 

n 

fn+l = X[kp 

2=1 

One conhrms that the fj satisfy the constraint (0)- (Note that the maximal sym¬ 
metry is evident in terms of the pfs since kj = pj — pn+i and kj — ki = pj — pi.) The 
proof that this is the correct form for the fj follows by constructing the Poisson 
brackets and collecting like exponentials. Conditions are quickly found that the fj 
must be particular products of differences between momenta. It is then seen that 
there are no additional requirements. 


(n -h 1) 
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Using and (|35|) in ( |^ and gives the explicit canonical transforma¬ 
tion between the open-chain n-body Toda Hamiltonian in the Chevalley basis 
and a Hamiltonian (^0|) which is independent of coordinates. The reduction to 
action-angle variables is essentially complete. From this point, one can attempt to 
construct a product of elementary canonical transformations which produces this 
full transformation. This has been done for the 3-body system and work is in 
progress on the n-body system. The value of such a product is that, when it is 
found in the quantum system, it allows the construction of integral representations 
of the eigenfunctions of the system. 
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